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Abstract 

We consider discrete nonlinear hyperbolic equations on quad-graphs, in particular 
on 1? . The fields are associated to the vertices and an equation Q(x±, X2, x$, X4) = 
relates four fields at one quad. Integrability of equations is understood as 3D- 
consistency. The latter is a possibility to consistently impose equations of the same 
type on all the faces of a three-dimensional cube. This allows to set these equations 
also on multidimensional lattices Z . We classify integrable equations with complex 
fields x, and Q affine-linear with respect to all arguments. The method is based on 
analysis of singular solutions. 



1 Introduction 



> 

The idea of consistency (or compatibility) is at the core of the theory of integrable systems. 
\q '. It appears already in the very definition of complete integrability of a Hamiltonian flow in 
the Liouville- Arnold sense, which says that the flow may be included into a complete family 
of commuting (compatible) Hamiltonian flows [1] . Similarly, it is a characteristic feature of 
soliton (integrable) partial differential equations that they appear not separately but are 
always organized in hierarchies of commuting (compatible) systems. The condition of the 
existence of a number of commuting systems may be taken as the basis of the symmetry 
approach which is used to single out integrable systems in some general classes and to 
classify them [18J. Another way of relating continuous and discrete systems, connected 
with the idea of compatibility, is based on the notion of Backlund transformations and 
the Bianchi permutability theorem [Hj. The latter developed into one of the fundamental 
principles of discrete differential geometry [T2] . 
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So, the consistency of discrete equations takes center stage in the integrability theater. 
We say that 

a d-dimensional discrete equation possesses the consistency property, if it may 
be imposed in a consistent way on all <i-dimensional sublattices of a (d + 1)- 
dimensional lattice 

(a more precise definition will be formulated below). As the above remarks show, the idea 
that this notion is closely related to integrability, is not new. In the case d = 1 it was 
used as a possible definition of integrability of maps in [23]. For d = 2 a decisive step 
was made in [10] and independently in [19]: it was shown that the integrability in the 
usual sense of soliton theory (as existence of a zero curvature representation) follows for 
two-dimensional systems from the three-dimensional consistency. So the latter property 
may be considered as a definition of integrability. It is a criterion which may be checked 
in a completely algorithmic manner starting with no more information than the equation 
itself. Moreover, in if this criterion gives a positive result, it delivers also the discrete zero 
curvature representation. 

Basic building blocks of systems on quad-graphs are quad-equations, which are equa- 
tions on quadrilaterals 

Q{x 1 ,x 2 ,x 3 ,x 4 ) = 0, (1) 

where the field variables Xi e CP 1 assi gned to the four vertices of the quadrilateral as 
shown in Figure CD On Z 2 equations of this type can be treated as discrete analogues 
of nonlinear hyperbolic equations. Boundary value problems of Goursat type for such 
systems were studied in [6]. 

Assumption. In this paper we assume that Q is affine- linear, i.e. a polynomial of 
degree one in each argument. This implies that equation ([T]) can be solved for any variable 
and the solution is a rational function of other three variables. 

X4 £3 



Q 



Xi x 2 

Figure 1: A quad-equation Q(x\, x 2 , x 3 , X4) = 0; the variables Xi are assigned to vertices 

The general idea of integrability as consistency in this case is shown in Figure EJ We 
put six quad-equations on the faces of a coordinate cube. The subscript j corresponds to 
the shift in the j-th coordinate direction. If one starts with arbitrary values x, Xi, x 2 , £3 
then the values Xi 2 , Xi 3 , x 23 are found from three equations on the left, front and bottom 
faces and the equations on the right, back and top faces yield, in general, three different 
values of x± 2 3. The system is called 3D-consistent, if these three values are identical for 
arbitrary initial data x, X\,x 2 , £3. 
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Figure 2: A 3D consistent system of quad-equations. The equations are associated to faces 
of the cube. 

In [3] we classified 3D-consistent systems of a particular type. The equations on all 
faces coincided up to the values of parameters associated with three types of the edges. 
Moreover, cubic symmetry was imposed as well as a certain additional condition called 
the tetrahedron property. A 3D-consistent system without the tetrahedron property was 
found in [13J. Later, this system was shown to be linearizable in (22]. In [2S] it was shown 
that the 3D-consistent equations classified in [3] satisfy the integrability test based on the 
notion of algebraic entropy. 

The consistency approach was generalized in various directions. Systems with the fields 
on edges lead to Yang-Baxter maps [23| [231 [20] ■ Quadrirational Yang-Baxter maps were 
classified in [3J. The 4D-consistency of discrete 3D-systems is related to the functional 
tetrahedron equation studied in [TTj, [16j [T5], [8]. 

In the present paper we classify 3D-consistent affine-linear quad-equations in a more 
general setup. The faces of the consistency cube can carry a priori different quad-equations. 
Neither symmetry nor the tetrahedron property are assumed. This leads to a general 
classification of integrable quad-equations. 

The outline of our approach is the following. 

a) By applying discriminant-like operators to successively eliminate variables one can 
descend from an affine-linear polynomial of four variables, associated to a quadrilateral, 
to quadratic polynomials of two variables, associated to its edges, and finally to quartic 
polynomials of one variable, associated to its vertices (Section [2]). 

b) By analysis of singular solutions, we prove that the biquadratic polynomials which 
come to an edge of the cube from two adjacent faces coincide up to a constant factor (see 
Section [3]). At this point an additional non-degeneracy assumption is needed. We assume 
that all the biquadratic polynomials do not have factors of the form x — c with constant 
c. (Examples of equations without this property are presented in Section [7]). 

c) This allows us to associate to each vertex of the cube a quartic polynomial in 
the respective variable; the admissible sets of polynomials are classified modulo Mobius 
transformations, each variable is transformed independently (Section HJ). 

d) Finally we reverse the procedure and reconstruct the biquadratic polynomials on 
the edges of the cube and the affine-linear equations themselves (Section [5]). 
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2 Affine-linear and biquadratic polynomials 

Our approach is based on the descent from the faces to the edges and from the edges 
to the vertices of the cube. In this Section we consider a single face and describe this 
descent irrespective of 3D-consistency. Let V™ denote the set of polynomials in n variables 
which are of degree m in each variable. We consider the following action of Mobius 
transformations on polynomials / G V™: 

M[f) {xi,...,x n ) = (c lXl + d.T . . . (c n x n + d n ) m f ( aiXi + h \ . . . , a " x " + M , 

\C\X\ + d\ c n x n + d n / 

where a^c/j — fejCj = Aj ^ 0. The operations 

T~*4 * T~*2 ~* ^1 ' $x,y(Q) = QxQy QQxyi $x(fo) = h x 2hh xx 

are covariant with respect to Mobius transformations. (The subscripts denote partial 
differentiation). More precisely: if Q G V}, h G Pf , then 

5 Xi>x .(M[Q\) = AAM^iQ)}, K{M[h]) = ^M[5 Xi {h)]. (2) 

Further on we will make an extensive use of relative invariants of polynomials under 
Mobius transformations. For quartic polynomials r G V\ these relative invariants are 
well known and can be defined as the coefficients of the Weierstrass normal form r = 
4x 3 — g 2 x — g%. For a given polynomial r(x) = r 4 x 4 + r 3 x 3 + r 2 x 2 + r%x + r they are given 
by (see e.g. 



g 2 (r,x) = ^-(2rr IV - 2r'r"' + (r") 2 ) = — (12r r 4 - 3rir 3 + r 2 ), 
4o LZ 

g 3 f r x ) = J—U2rr"r IV - 9(r') 2 r IV - 6r(r'") 2 + 6rVV" - 2(r") 3 ) 
3456 

(72ror2r 4 — 27r 2 r 4 + 9rir 2 r^ — 21r^r2^ — 2r 2 ). 



432 

Under the Mobius change of x — x\ these quantities are just multiplied by the constant 
factors: 

g k (M[r],x) = Al k g k (r,x), k = 2, 3. 
For biquadratic polynomials h G V 2 , 

h(x, y) = h 22 x 2 y 2 + h 21 x 2 y + h 20 x 2 + h 12 xy 2 + h u xy + h w x + h 02 y 2 + h m y + h 00 , (3) 

the relative invariants are 

t 2 (hj x, y^j 2hh XX yy 2h x h X yy 2hyh XX y ~\- 2h xx hyy ~\- fa x 



xy 

11) 

" h x h xx \ I h 22 h 2 \ h 2 Q 



Sh 00 h 22 - Ah 01 h 21 - 4/ii /ii2 + %h 02 h 2Q + h 2 



i 3 (h,x,y) = -det I h y h xy h xxy I = det I h X2 h u h w 
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■ hyy h X yy fl XX yy J \ ^02 ^01 ^"00, 
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Notice that i 3 can be denned as well by the formula 

-4i 3 (h,x,y) = 5 X) y(5 X) y(h))/h. 

Under the Mobius change of x — X\ and y = x 2 , 

i k (M[h],x, y) = A*A*i k (h, x,y), k = 2, 3. 

The following properties of the operations 8 xy , 5 X are proved straightforwardly. 

Lemma 1. For any affine-linear polynomial £2, £3, £4) G V\ there holds: 

s x 3 (8 xl 2 (Q)) = S X2 (5 xliX3 (Q)), 
ik(S XljX2 (Q),x 3 ,x i ) =i k (8 X3tX4 (Q),xi,x 2 ), k = 2,3. 

For any biquadratic polynomial h(xi,X2) G V\ there holds: 

gk(S Xl (h),x 2 ) = g k (S X2 (h),Xx), k = 2,3. 



(4) 
(5) 



(6) 



Denote = h? % = S XktXl (Q) where k,l} = {1,2,3,4}. Then Lemma [H implies the 
commutativity of the diagram 



r 4 (x 4 ) 



S X1 



J X4 



n(xi) 



h u (x 3 ,x 4 ) 



J X4 



7 II, °X2,X3 r\l \ °X1,X4 ,93/ \ 

h^[x u x A ) i Q{x 1 ,x 2 ,x 3 ,x 4 ) > h za (x 2 ,x 3 ) 



J X3,X4 



r 3 {x 3 ) 



'£■3 



---- h 12 (x!, x 2 ) > r 2 (x 2 ) 



(7) 



Moreover, biquadratic polynomials on the opposite edges have the same invariants i 2 ,i 3 , 
and all four quartic polynomials r, have the same invariants g 2 ,g 3 . This diagram can be 
completed by the polynomials h 13 , h 2A corresponding to the diagonals (so that the graph of 
the tetrahedron appears), but we will not need them. The introduced polynomials satisfy 
also a number of other interesting relations. 

Lemma 2. For any affine-linear polynomial Q{x%, x 2 , x 3 , X4) G V\ and with the notations 
h^(xi,Xj) = o~x k ,xi{Q) G V 2 , the following identities hold: 

ti 3 (h 12 ,x u x 2 )h u = det I fcg hf lX2 £ X2 



h 12 h 12 £ 

IL X2X2 IL X 1 X2X2 ^^2X2 



when I h%..lr> - hfjr" • HPh^; 



X 3 '"X 3 



x 3 x 3 > 
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Q Qx± Qx3 

h u h u _ ^23 = P g ; P = det | Q X2 Q X1X2 Q X2X3 | g V\ ; (9) 

QxiX4 QxgX4 



2_Qan h x 2 h 3A -h x \h 23 + h 23 h 3A -h 2 x \h 3A 

Q ft 12 ft 34 - h 14 h 23 



(10) 



Identity ([8]) shows that ft 14 can be expressed through three other biquadratic poly- 
nomials (provided ^(ft 12 ) ^ 0). Identity ([9]) defines Q as one of the factors in a simple 
expression built from ft*- 7 . Finally, differentiating fflUj) with respect to x 2 or ^4 leads to a 
relation of the form Q 2 = F[h 12 , ft 23 , ft 34 , ft, 14 ], where F is a rational expression on ft 4 - 7 and 
their derivatives. Therefore, if the biquadratic polynomials on three edges (out of four) is 
known, then Q can be found explicitly. Of course it is seen from Lemma [2] that not any 
set of three biquadratic polynomials comes as ft*- 7 from some Q &V\. 

Biquadratic polynomials ft*- 7 for a given Q G V\ are closely related to singular solutions 
of the affine-linear equation 

Q{x 1 ,x 2 ,x 3 ,x 4 ) = 0. (11) 

The polynomial Q G V\ is assumed irreducible (in particular, Q Xi ^ 0: otherwise the 
polynomial Q should be considered as reducible, since under the change of the variable 
Xi 1— > 1/xi it turns into x^Q). Obviously, equation ( TTTT) can be solved with respect to any 
variable: let Q = p(xj,Xk,Xi)xi + q(xj,Xk,Xi) then = —q/p for the generic values of 
Xj, Xk, X[. However, Xi is not determined if the point (xj,Xk, x{) lies on the curve in (CP 1 ) 3 

Si\ p(xj,x k ,xi) = q(xj,x k ,Xi) = 0, Q = pxi + q. (12) 

The projection of this curve onto the coordinate plane (j, k) is exactly the biquadratic 
h? k = pq Xl -p Xl q = 0. 

Definition 1. A solution (xi, x 2 , x$, X4) of equation (TTTT) is called singular with respect 
to Xi , if it satisfies also the equation Q x .(xi,x 2 , x 3 , £4) = 0. The curve is called the 
singular curve for Xj. 

Lemma 3. If a solution (xi, x 2 ,x^, X4) of equation < f77]j zs singular with respect to X{, then 
W k = W l = h kl = 0. Conversely, ifW k = for some solution, then this solution is singular 
with respect to either Xi or x\ . 

Proof. Since ft jfc = Q Xi Q Xl — QQx,,x v the equations ft jfc = and Q Xi Qx t = are equivalent 
for the solutions of equation Q = 0. □ 

We will use the following notion of non-degeneracy for biquadratic polynomials. 

Definition 2. A biquadratic polynomial h(x,y) G V 2 is called non- degenerate, if no 
polynomial in its equivalence class with respect to Mobius transformations is divisible by 
a factor x — c or y — c (with c = const). 
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Figure 3: A 3D consistent system of quad-equations. The equations are associated to faces 
of the cube: A and A on the bottom and on the top ones, B and B on the front and on 
the back ones, C and C on the left and on the right ones, respectively. 

According to this definition, a non-degenerate polynomial h(x,y) G V\ is either irre- 
ducible, or of the form {a.\xy + fox + 71?/ + Si) (a 2 xy + fox + r ) 2 y + #2) with o^Si 7^ foji. In 
both cases the equation h = defines y as a two-valued function on x and vice versa. On 
the other hand, for example the polynomial h(x, y) = x — y 2 (considered as element of V 2 ) 
is, according to Definition [2], a degenerate biquadratic, since under the inversion x 1— > 1/x 
it turns into x(l — xy 2 ). 

A fundamental role in our studies will be played by the following notion. 

Definition 3. An affine-linear function Q G is said to be of type Q, if all four its 
accompanying biquadratics h? k e V\ are non-degenerate, and is said to be of type H 
otherwise. 



3 3D-consistency and biquadratic curves 

Consider the system of equations 

A(x, xi,x 2 , x u ) = 0, A(x 3 , xi 3 , x 2 3, X123) = 0, 

B(x, xi,x 3 , x i3 ) = 0, B(x 2 , X12, x 2 3, X123) = 0, (13) 

C(x, x 2 , x 3 , x 23 ) = 0, C(xi, X12, xi 3 , X123) = 

on a cube, see Figure El The functions A,...,C are afline linear (from V\) and are 
not a priori supposed to be related to each other in any way. We will use the notation 
A % i = S XkjXl A for the accompanying biquadratic polynomials. 

Theorem 1. Let all six functions A,...,C be of the type Q, and let equations 173]) be 
3D- consistent. Then 
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1 ) for any edge of the cube, the two biquadratic polynomials corresponding to this edge 
( coming from the two faces sharing this edge ) coincide up to a constant factor; 

2) the product of these factors around any vertex is equal to —1, for example 

^0,1^0,3^0,2 + ^0,2^0,1^0,3 = Q . (14) 

3) the system 173]) possesses the tetrahedron property dx^/dx = 0. 
Proof. Elimination of X12, X13 and X23 leads to equations 

F{x, x±, X2, X3, X123) = A xl3jX23 BC — A X23 B Xl3 C — A Xl3 BC X23 + AB Xl3 C X23 = 0, 
G(x, Xi, x 2 , X3, X123) = B xl2iX23 AC — B X23 A Xl2 C — B Xl2 AC X23 + BA Xl2 C X23 = 0, 
H(x, Xi, x 2l X3, X123) = C Xl2jX13 AB — C X13 A X12 B — C Xl2 AB Xl3 + C A X12 B X13 = 0. 

Here the numbers over the arguments of F, G, H indicate the degrees of the right hand side 
in the variables (the degree is understood in the projective sense, as in the example at the 
end of the previous Section). Due to 3D-consistency, the expressions for X123 as functions 
of x, X\ } X2, X3 found from these equations, coincide. Therefore these factorizations hold: 

F = f(x,x 3 )K, G = g(x,x 2 )K, H = h{x,x 1 )K, K = K(x, x Xl x 2 , x 3 , x 123 ), (15) 

where /, g, h are some polynomials of degree 2 in the second argument. The degree in x 
remains to be determined. 

Let the initial data x,xi,x 2 be free variables, and X3 chosen to satisfy the equation 
f(x, X3) = 0. Then F = 0, and thus the system B = C = A = does not determine 
the value of Xi 2 3- Moreover, the equation B = can be solved with respect to x i3 since 
otherwise the initial data must be constrained by equation B 0,1 (x, xi) = 0. Analogously, 
the equation C = is solvable with respect to £23. Therefore, the uncertainty appears from 
the singularity of equation A = with respect to £123. Hence, the relation A 3 ' 13 (x 3 , £13) = 
is valid. In virtue of the assumption of the theorem, x 13 is a (two- valued) function of x 3 
and does not depend on x\. This means that the equation B = is singular with respect 
to x% and therefore B 0,3 (x, x 3 ) = 0. Analogously, C 0,3 (x, x 3 ) = 0. 

Thus, we have proven that if X3 = <p(x) is a zero of the polynomial / then it is also a 
zero of the polynomials B 0,3 , C ' 3 . If one of these three polynomials is irreducible, then 
this already implies that they coincide up to a constant factor. If the polynomials are 
reducible this may be wrong since it is possible that f = a 2 , B ' 3 = ab, C 0,3 = ac, where 
a,b,c are affine-linear in x, X3. In any case we have deg^, / = 2 and this is sufficient to 
complete the proof. 

Indeed, this implies deg^ K = 0, that is the tetrahedron property is valid. In turn, this 
implies the relation ( TUP) , as was shown in [3]. Recall this calculation: let us rewrite the 
system fll3p in the form 

x 12 = a (x, xi, x 2 ), x 13 = b(x, xi, x 3 ), x 23 = c(x, x 2 , x 3 ), 
xu3 = d(xi, x 2 , x 3 ) = a(x 3 , x 13 , x 23 ) = b(x 2 , x 12 , x 23 ) = c(xi, x 12 , x 13 ) 
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and find by differentiation 



d Xl 


^xi3 b X i i 


d X2 


= CLx 23 Cx2 ) 





^xi^bx 




dx\ 




dx3 







bxi2^x 


-\- b X2S c X j 


dx2 


= C Xl2 CL X2 , 


dx% 


— c Xl3 bx 3 i 





Cx\ 2 ^x 


~\~ c Xl3 b x - 



These equations readily imply the relation 

Q j x2bx\Cx'i ^x\bx 3 c X2 0. 

The latter is equivalent to by virtue of the identity a X2 /a Xl = A 0,1 /A 0,2 . The variables 
in equation separate: B°> 3 /C°> 3 = -A°' 2 /C°> 2 ■ B^/A * 1 , so that B°' 3 /C ' 3 may only 
depend on x. In view of the assumption of the theorem this ratio is constant. □ 

There exist 3D-consistent systems whose equations are not of type Q. The statements 
of Theorem [I] may or may not be valid for such a system, as the following examples 
demonstrate. 

Example 1. The simplest 3D-consistent equation is the linear one 

x I x 2 I x j I x £ j — . 

In this case, all the biquadratic polynomials are equal to 1, so that the statement 1) is 
fulfilled and statement 2) is not. Since 2) is a consequence of the tetrahedron property 3), 
the latter cannot hold either. Indeed, 

X123 = 2x + X 1 + X 2 + X 3 . 

The factor / in this example is also equal to 1, but this is a coincidence, destroyed by 
Mobius changes of variables. Indeed, in this case deg x K = 1, and after the inversion of 
all variables xj — > 1/xj we come to / = xx 2 , while B 0,3 turns into x 2 x\. 

Example 2. The Hietarinta equation [13] 

{x-e U) ) ( Xij - o® ) ( Xi - o (i) ) ( Xj - e (i) ) - (x - e (i) ) ( Xij - o (i) ) (z ( - ) ( Xj - o {j) ) = (16) 
is 3D-consistent, but the statement 1) does not hold: 

£0,3 = (e (3) _ (l) )(o (l) _ (3) )(x _ e (3) )(x _ e (l) )(x3 _ e (l) )(x3 _ o(3))) 
C 0,3 = (e (3) _ (2) )(o (2) _ (3) )(x _ e (3) )(x _ e (2) )(x3 _ ^2)^ _ ( 3)) _ 

The factor / is proportional to (x — e^ 3 ^)(a;3 — eW)(i 3 — e^). Consequently, deg x K = 1 
and the tetrahedron property does not hold. 

Example 3. Probably the best known example of a 3D-consistent system is given by the 
discrete potential KdV equation 

(x - Xij)(xi - xj) + a {i) - a {j) = 0. (17) 

In this case all statements of the theorem are valid, in spite of the degeneracy of biquadrat- 

ICS" 

B°> 3 = - a®, C ' 3 = a® - a®, f = 1. 

(Recall that the degree is understood in the projective sense. Under the inversion these 
polynomials turn into x 2 x 2 .) 
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Example 4. Equation fl^D 

Q(x, xi, x 2 , x 12 ; a (2) ; 5) 

= a^\x — x 2 )(xi — xi 2 ) — ot^(a; — Xi 

is consistent not only with its own copies (see [3] and Theorem H] below) , but also with 
the linear equations. Namely, the system formed of the equations 



Q(x, xi, X12, x 2 , a (2) ; 5) = 0, 



X13 - x 3 = xt - x, x 23 - x 3 = x 2 - x 



and their copies on the opposite faces, is 3D-consistent. In this case the edge (x, x 3 ) carries 
the polynomials 

£0,3 = ^0,3 = _ 1; f = L 

However, in contrast to the previous example, the tetrahedron property is not valid and 
deg x K = 2. This means that the polynomial / is not biquadratic and its image under 
inversion is x 2 . Moreover, the biquadratic polynomials corresponding to the edge (x,xx) 
do not coincide: 

A > 1 = aP\ a V>-a®)({xi-x) 2 -6( a M) 2 ), B ' 1 = -1, h = 1. 

We see in this example that it is possible that some of the biquadratic polynomials satisfy 
the assumption of the theorem and the others do not. 



4 Classification of biquadratic polynomials 

Diagram (J7|) suggests an algorithm for the classification of affine-linear equations Q = 
modulo Mobius transformations. The first step is to use Mobius transformations to bring 
the polynomials r^(xj) associated to the vertices of the quadrilateral into canonical form. 
According to formulas ([2]), 

5 Xl (S Xj , Xk (M[Q])) = A 2 A 2 A 2 M[M^, Xfc (Q))] = ^M[r<], 

i 

where C = Af A 2 A|A|. Since the polynomial Q is defined up to an arbitrary factor, we 
may assume that Mobius changes of variables in the equation Q = induce transformations 

i 

of the polynomials tv This allows us to bring each into one of the following six forms: 

r = (x 2 — l)(A; 2 x 2 — 1), r = x 2 — 1, r = x 2 , r = x, r = 1, r = 0, 

according to the six possibilities for the root distribution of r: four simple roots, two 
simple roots and one double, two pairs of double roots, one simple root and one triple, one 
quadruple root, or, finally, r vanishes identically. Note that in the first canonical form it 
is always assumed that k 7^ 0, ±1, so that the second and third forms are not considered 
as particular cases of the first one. 
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Not every pair of such polynomials is admissible as a pair of polynomials at two adja- 
cent vertices, since the relative invariants of the polynomials of such a pair must coincide 
according to We identify all admissible pairs, and then solve the problem of recon- 
struction of the biquadratic polynomial (E]) by the pair of its discriminants 

h 2 — 2/iZi,,,, = 7*1(2;), S x (h) = h 2 — 2hh xx 



5y{h) 



J yV"J '"y ^<«,Uyy 

which is equivalent to a system of 10 (bilinear) equations for 9 unknown coefficients of the 
polynomial h. 

Theorem 2. Biquadratic polynomials with the pair of discriminants (7*1(2;), 7*2(7/)) in 
canonical form exist for the following pairs, up to the permutation of x, y: 







y 2 - 1 y 2 


y 1 





(x 2 - l)(k 2 x 2 - 
x 2 -l 


-1) 


+ 












+ + 






x 2 






+ 






X 








+ + 




1 








+ 


+ 













+ 



These polynomials h and their relative invariants i 2 ,i3 are: 
(r(x),r(y)), r(x) = (x 2 — l){k 2 x 2 — 1) : 

2Axy — x 2 
3(k 2 a 2 + a" 2 ) - k 2 - 1, 4z 3 

,2n 



h 

1 2 

2 x „;2 



1 i 1 2 2 2 2 

— ( k ex. x y 



2 1 2 
y +a 



A 2 = r(a), 
A(k 2 a - a" 3 ); 



x 2 - 5, y 2 - 5) : h 
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1-a 2 
1 + 10a 2 + a 4 

«3 



;i - a 2 ) 2 
1 



[x< + y< 
a 2 (l + a 1 



'l-a 2 ) 3 



1 + a 2 5(1 
-xy H 



a 



1-a 2 



4a 



(19) 



(20) 



(x - y) A 



a .a 
-{x + y) + — 



1 2 



«3 



1 



(x,y):h-^ „, 2 ^, y/ , 4 , ^- 4a2 , -o 32a 3 

(x 2 , y 2 ) : h = \x 2 + fixy + uy 2 , ji 2 — 4Az/ = 1, i 2 = 1 + 12Az/, i 3 = —X\iv\ 
Xx 2 y 2 + fixy + u, /1 2 — 4Az/ = 1, i 2 = 1 + 12Azv, i 3 = A/iz/; 



h 



(1,1): /i = A(x ± y) 2 + /i(x ± y) + i/, 
(0, 0) : h = (xxy + Ax + jj,y + z/) 



(x 2 -l,y 2 ) 



h = ay ±xy + — , 



/i — 4Az/ = 1, %2 
i 2 = 12(xz/ — A/i) 2 , 

i 2 = 1, i 3 = 0; 



>3 



12A 2 , 23 = T-2A 3 ; 
= 2(xz/ - A/i) 3 ; 



(x, 1) : /i = ±-{y -a) 2 Tx 



H = 0, i 3 = 0; 
1,0): /i = Xy 2 + fiy + z/, /i 2 — 4Az/ = 1, 2 2 = 0, i 3 



0. 



(21) 

(22) 
(23) 
(24) 
(25) 

(26) 

(27) 
(28) 



Proof. The list is obtained by solving the system ( Tl8l) for different canonical pairs (7*1,7*2). 
The exhaustion of cases is shortened if we notice that g 2 7^ 27gf in one case only, and that 
the relative invariants for the polynomial 7*1 = ax 2 + bx + c are 12g 2 = o 2 , 216(73 — — o 3 , 



so that the second polynomial must be of the form r 2 
pair (x, 0) turns out to be empty. 



ay +by + c. The solution for the 

□ 
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5 Classification of affine-linear equations of type Q 



It is important to note that after bringing the polynomials rj(xj) into canonical forms, one 
still has some freedom. Namely, one can use Mobius transformations which do not change 
the form of r to further simplify the biquadratics h and the affine-linear equation Q. In 
particular, the list of Theorem [2] is slightly more detailed than the list of biquadratics 
modulo Mobius transformations. 

Indeed, the polynomial (I2"2"j) turns into (1231) under the inversion of x; the change x h- > —x 
allows to fix the signs in the polynomials (I24p . ( 1261) ; in the case (1271) . the sign is fixed by 
the change x ^ —x, y ^ iy; the polynomials (1251) . (1281) admit a further simplification. 

However, a transformation of any one of the four variables for a quadrilateral influences 
biquadratic polynomials on two edges adjacent to the correspondent vertex, and therefore 
it cannot a priori be guaranteed that all four biquadratics can be brought to some definite 
form simultaneously. For example, if each vertex corresponds to the polynomial = x 2 
then the edges may correspond to the polynomials either of the form ( 1221) or of the form 
(123]) . We do not know a priori that these polynomials can be always brought into the 
same form (even with different coefficients). Actually, this is possible, as the proof of the 
following theorem shows. 

The next step is the reconstruction of the affine-linear polynomials from the biquadratic 
ones. Since our goal is only the classification of systems of type Q equations, we will not 
solve this problem in its full generality. We leave aside the cases ( )26l) . ( |271) and ( 1281) . 
since the corresponding biquadratics are degenerate. For the same reason, the additional 
restrictions on the values of parameters are imposed: \v ^ in the cases (I2"2"j) . (I2"3"j) . A 7^ 
in the case (12^1) . and kv — A/i 7^ in the case ( 125*1) . 

Theorem 3. Any affine-linear equation of type Q is equivalent, up to Mobius transforma- 
tions, to one of the equations in the following list: 

sn(a) sn(/3) sn(a + (3)(k 2 xix 2 x 3 x 4 + 1) — sn(a)(x\x 2 + ^3^4) 

— sn(/3)(xix 4 + x 2 x 3 ) + sn(a + /3)(xiX 3 + x 2 x 4 ) = 0, (29) 
(a - a~ 1 )(x 1 x 2 + x 3 x 4 ) + ((3 - (3~ 1 )(x 1 x 4 + x 2 x 3 ) - (a/3 - a~ 1 /3~ 1 )(x 1 x 3 + x 2 x 4 ) 

+ -(a - a~ v ){f3 - /? _1 )M - a" 1 ^" 1 ) = 0, (30) 

a(xi - x 4 )(x 2 - x 3 ) + (3{xi - x 2 )(x 4 - x 3 ) 

- a/3(a + (3)(x! + x 2 + x 3 + x 4 ) + a(3(a + (3)(a 2 + a/3 + (3 2 ) = 0, (31) 
a(xi — x 4 )(x 2 — x 3 ) + (3(xi — x 2 )(x 4 — x 3 ) — 5a(3(a + (3) = 0. (32) 

Proof. Let the polynomials h 12 , h 23 , h 3A and h lA be of the form ( fl9l) . with the parameters 
(a, A), (f3, B), (a, A) and ((3,B), respectively, lying on the elliptic curve A 2 = r(a). The 
relative invariants i 2 , i 3 of h 12 and h 3A must coincide as a corollary of (jSJ), and it is easy to 
check that this condition allows only the following possible values for (a, A): 

(a, A), (-a, -A), ^(a,-A), —(-a, A) 

and analogously for (/3, B). At first glance, it seems like we had to examine 16 quadruples 
of h % i , but actually the situation is much more favorable. Indeed, according to (j2J), a 
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Mobius change of variables in the equation Q = yields 

5 XkjXl (M[Q}) = A k \M[S ahm (Q)] = ^M[h% 

where C = A1A2A3A4. Since Q is only defined up to a multiplicative constant, we may 
assume that a Mobius change of variables induces transformations 

h ij ^ ——M[h ij ] 
A., A,- L J 

on the biquadratic polynomials h tJ . In particular, if 

1 A 



h 34 = h(x 3 , x 4 , — a, — A) or h 34 = h(x 3 , x 4 , 



-hi— x 3 , x 4 ; — a, — A), resp. — kxlh( - — ,x 4 ; 

\kx-i 



ka ka 2 . 

then the corresponding one of the Mobius transformations x 3 1— > — a; 3 or x 3 1— ► 1/ (fea^) will 
change /i 34 to 

1 1 _ A 

both of which coincide with h(x 3 , x 4 , a, A) due to symmetries of the polynomial (Tl9l) . Thus, 
performing a suitable Mobius transformation of the variable 0:3 (which does not affect the 
polynomial r(x 3 )), we may assume without loss of generality that (a, A) = (a, A). After 
that, the polynomial h 14 is uniquely found according to formula (jHJ), and it turns out that 
the equality ($, B) = B) is fulfilled automatically. Thus, the change of one variable 
allows to achieve the equality of the parameters corresponding to the opposite edges of 
the square. A direct computation using formula ( JTU1) yields the equation 

a(5 r ){k 2 XxX2X^x i + 1) + a(xix 2 + x 3 x 4 ) + fi{x\X± + x 2 x 3 ) + 7(^1X3 + x 2 x 4 ) = 0, 

where 7 = (aB + /3A)/(k 2 a 2 /3 2 — 1), and the final change a — > sn(a), A — > sn'(a) and 
analogously for f3 brings it to the form (1291) . 

Also in the other cases, suitable Mobius changes of the variables x 2 , x 3 , x 4 allow 
us to bring the polynomials into the form h 12 = h(xi,X2,a), h 23 = h(x2, X3, (3), h 3A = 
h(x 3 ,X4,a). Moreover, a direct computation using formula (jHJ) proves that also h 14 = 
h(x±, X4, (5). Then the answer is found by use of ( ITOl) . 

To give a few more details, polynomials ( 120|) give rise to equation ( |30|) . In this case 
equations (jSJ) imply that the parameters a of the polynomials h 12 and h 34 differ at most by 
sign. This is compensated by the change x 3 — > —x 3 which is possible due to the symmetry 
h(x,y,a) = —h(—x,y, —a). 

In the cases ( )22l) . ( J23l) . the appropriate scalings and, if necessary, inversions of the 
variables 22, £3, x 4 allow us to bring h 12 , h 23 , h 34 into the form (1201) without the constant 
term; therefore we simply obtain this case at 5 = 0. 

Polynomial (j2T|) corresponds to equation flHTl) . This is the most simple case since the 
parameters are fixed already by condition (j3J). 

In the case (IMj) . appropriate shifts and, if necessary, changes of sign of the variables 
£2, 2^3, x 4 allow us to bring h 12 , h 23 , h 34 into the form 2h(x,y,a) = a^ 1 (x — y) 2 — 5a 
with 5 = 1. Analogously, in the case ( 1251) appropriate Mobius transforms of the general 
form bring h 12 , h 23 , h 34 into the same form with 5 = 0. In both cases, the invariants are 
%2 = 3a~ 2 , 4i 3 = a~ 3 , therefore the parameters of h 12 and h 34 coincide and no further 
changes are necessary. The resulting equation is (132]) . □ 
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6 Classification of 3D-consistent systems of type Q 



Theorem [T] provides very strong necessary conditions for 3D-consistency in the case when 
all equations are of type Q. This will allow us to classify such systems in this section. In 
this final step we have to arrange the obtained equation around the cube and to choose the 
parameters in such a way that the condition (1141) is fulfilled. The effect of this condition 
may be a change of sign or an inversion of one of the parameters. 

In the following Theorem we return to the notation of the variables and parameters 
according to the shifts on the lattice. The ordering of the equations corresponds to the 
previous Theorem, and we name these equations as in [3]. 

Theorem 4. Any 3D-consistent system (Q3|) of type Q is, up to Mobius transformations, 
one of the following list: 

sn(a^) sn(a ( -^) sn(a^ — a^)(k 2 xXiXjXij + 1) + sn(a^)(xxi + XjXij) 

— sn(a^)(xxj + XiXij) — sn(a® — a^)(xxij + XiXj) = 0, (Q^) 

1 \ , s /«W a ij) ' 



\( \ ( (J) \( \ ( oi yj <\ 



6/ ,« 1 W 1 \ fa® a® 



. a< " ~ r w - soy J - |» ) = °- «w 

Oi ^X X ' j ^ (^X >i X ij J 01 (^X X 'i J (^X j X J I 01 Gt ^ 01 01 ^ (^X \ X^ \ X j \ X^j ^ 

- a®a®(a® - a^)((a») 2 - a {i) a {j) + (a ij) ) 2 ) = 0, (Q 2 ) 
- Xj)(xi - x,ij) - a (jf) (x - Xi)(xj - x^) + 5a^a^\a® - a (jf) ) = 0. (Qi) 

Proof. First of all, note that equations of the different types (I29~ i) -( l32"l) cannot be consistent 
with each other since the corresponding singular curves are different. In particular, the 
parameters k 2 in the case ( 1291) and S in the cases ( 130|) . ( |32j) must be the same on each face 
of the cube. Moreover, each equation of the list possesses the square symmetry, that is, it 
is invariant with respect to the changes (xi <-> x 2 , ^3 *-* x i) an d O^i ^ x 3 , a ^ (3). 

Therefore, the equations on all faces may differ only by the values of a and f3. Consider 
the equations corresponding to three faces meeting in one vertex, say x: 

Q(x, xx,x 2 , X12, a, P) = 0, Q(x, x 2 , x 3 , x 23 , P, 7) = 0, Q(x, x 3 , x x , x 13 , 7, a) = 0. 

Let 

S X2)Xl2 Q{x, xi, x 2 , xu, a, (3) = K(a,j3)h(x,xi,a). 
Then, due to the symmetry, 

S Xl , Xl2 Q(x, xx,x 2 , x 12 , a, 0) = K0, a)h(x, x 2 , (3) 

and, according to the Theorem [TJ the parameters must be related as follows: 

h(x,x l7 a) h(x,x 2 ,(3) ~ h(x,x 3 ,i) 

— — =m(a,a), — ~- = m{p,p), — — = 7727,7, 

h{x,x 1 ,a) h(x,x 2 ,(3) h{x,x 3 ,j) 

K(a,/9)K(/3,7)K(7,fi) 



«(^,«)«(7»^)«(«>7) 



m(a, a)m(f3, /3) 771(7, 7) = — 1. 
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In the case ( 1291) . a direct computation proves that hi(a, (3) = 2 sn(a) sn(/3) sn(a + (3) and 



h(x,y,a) 



1 



2 sn(aO 



(fc sn (a)x y + 2sn'(a)xy — x —y + sn (a)), 



therefore a may take the values ±a and analogously for /3, 7. Obviously, up to a change 
of the numeration, two cases are possible: 



a 



-a, (3 = -0, 7 = -7 



or 



a = a, (3 = j3, 7 = —7. 



Moreover, this is actually only one case since we can make the change (a,/?) — > (—a, —p) 
which keeps the equation Q(x, x±, x 2 , Xi 2 , a, $) = invariant, as one can easily see from 
( 1291) . It is not difficult to check that we can always adjust the signs on the whole cube as 



in the system QQ 4 p . 

Consider now the case (I30I) . Here 



(3) 

h(x,y,a) = - 



(l-a 2 /5 2 )(l-« 2 )(l-/3 2 ) 



a 



a' 



(x 2 + y 2 ) 



a 2 [3 2 
1 + a 2 



Q- 



-xy + 



a 2 )S 



4a 



and a = a or a = l/a. Taking into account the invariance of equation (l30i) with respect 
to the simultaneous inversion of a, (3, we can set, without loss of generality, 

5 = l/a, P = l/{3, 7 = 1/7 



resulting in the system QQ 3 [ ). In the cases ( 1311) . (1321) we have respectively 

/c(a,/3) = -4af3(a + (3), h(x,y,a) = —(x - y) - ~(x + y) + —, 



1 2 a6 

— (x — y) 

2a v yj 2 



n(a, (3) = — 2a(3{a + (3), h(x,y,a) 
and we may set a = —a, [3 = —j3, 7 = —7 exactly as before. This leads to the systems 

ra>, ra- □ 



The master equation <\Q^ of the list was first derived in [2] and further investigated in 
[5J. A Lax representation for (Q2) was found in [T5] with the help of the method based on 
the three-dimensional consistency. Equations QQjJ and dQ 3 | ^ =0 ) go back to [21]. Equations 
(^2) and (^3 5 = i ) first appeared explicitly in [3]. 



7 Examples of type H systems 

In contrast to the type Q systems, the systems of type H can be considered as "degener- 
ate" . Their classification seems to be a rather tedious task. Presently we cannot suggest 
any effective procedure to solve this problem. On the other hand, the examples given in 
the Section [3J demonstrate that this class should not be just neglected as "pathological". 
Indeed, the discrete KdV example (ITT]) suggests that in some cases the degeneracy of the 
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biquadratics is just an unessential coincidence which does not spoil the integrability prop- 
erties of an equation. Here we consider several more examples of this kind, corresponding 
to the cases (122]). (1231) at A// = 0. (123j) at A = 0. and (1251) at kv - A// = which were 
excluded in the previous Section. It turns out that if we apply the same algorithm in these 
cases (in spite of the fact that there is no justification for this) then the list H from our 
previous paper [3J will be reproduced: 

a {i) (xxi + xjXij) - a u) (xxj + XiXy) + 5((a {i) ) 2 - (a (i) ) 2 ) = 0, (H 3 ) 
(x - x^ixi - Xj) + (a U) - ct (i) )(x + Xi + xj + xtj) + (a ij) ) 2 - (a (i) ) 2 = 0, (H 2 ) 
(x - x^ixi - Xj) + a u) - a {i) = 0. (#i) 

One may check directly that all statements of the Theorem [T] remain valid for these equa- 
tions, in spite of the degeneracy of the biquadratics. 

Considering the asymmetric cases ( 1261) . ( 1271) . ( 1251) with the different polynomials asso- 
ciated to the different vertices, one finds that the following variants are possible, up to the 
permutations. (There is clearly no distinction between edges and diagonals when we are 
dealing with a single equation.) 

/ 2 i 2 2 2\ I 2 i 2 i 2 2\ (2 i 2 i 2 i 2\ 

[X^ 1, X 2 , X 3 , X^), [X^ *-j% 2 J -' -^3' y^l *-i x 2 ^-i^Z *-i x A)i 

(xi, 1,1,1), (xi,x 2 , 1, 1), (xi,x 2 ,x 3 , 1), 

(1,0,0,0), (1,1,0,0), (1,1,1,0). 

A direct check shows that the variants of type ( r2 ^" r4 | r i( x 3) \ arg rea ^ za ]3j e anc [ j eac [ ^ Q 



the following list of 3D-consistent equations: 



ri(xi) r 2 (x 2 ) 



££2X4 \ 



a(xix 2 + X3X4) - fl{x\x± + x 2 x 3 ) + (a 2 — /3 2 )[5 H — I 

V ap / 

(xi - x 3 )(x 2 - x 4 ) + (/3 — o;)(a;i + x 2 + x 3 + x 4 ) + (3 2 - a 2 

+ e(/3 - a)(2x 2 + a + f3)(2x A + a + (3) + e{j3 - af = 0, 
(x\ - xs)(x 2 - x 4 ) + ((3 - a)(l + £X2X 4 ) = 0. 



(#f) 



This list can be considered as a deformation of the list H which corresponds to the case 
e — 0. However, we use the notation with cyclic indices rather than shifts since due to 
the lack of symmetry the arrangement of the equations on the faces of a cube requires a 
more explicit description (see below). Note that in (7Tf) the polynomial 1 + ex 2 x± can 
be replaced by the polynomial kx 2 x a + (J,(x 2 + x 4 ) + v with arbitrary coefficients. The 
corresponding biquadratic polynomials and their discriminants are given in the following 
table (up to multiplication by a suitable constant, Q — * n(a, (3)Q): 



(HI) 



h(xi,x 2 ) 



ri(xi) r 2 (x 2 ) 



x\x 2 + ea x x\ + 5a x\ — 45e 
X\ + x 2 + a + 2e(x 2 + a) 2 1 — 8exi 
1 + ex\ -As 



y. 2 

'< 2 
1 
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Each of these equations possesses the rhombic symmetry 

Q(xi, x 2 , x 3 , x 4 , a, (3) = -Q(x 3 , x 2 , xi, x 4 , (3, a) = -Q{x x ,x^ x 3 , x 2 , (3, a), 

but not the square symmetry since the vertices Xi,x 2 correspond to polynomials with 
zeroes of different multiplicities. The equation is 3D-consistent on the black-white lattice 
i + j + k (mod 2). That is, each face must carry a copy of the equation in such way that 
the parameters on opposite edges coincide and the vertices x, x± 2 , Xi 3 , X23 are of the same 
type (here we switch again to the notation where indices denote shifts, as in Figure [3]): 

Q(x,x i ,x ij ,x j ,a { - t \a {: ' ) ) = 0, Q{x ih , x k , x jk , x 123 , Qt w , = 0, {i,j, k} = {1,2,3}. 

Obviously, the equations on opposite faces of the cube do not coincide, but the equation 
may nevertheless be extended to the whole lattice Z 3 . The tetrahedron property is fulfilled. 

Finally, we notice that it is also possible to combine equations with the square and 
trapezoidal symmetry. Consider equation (Qi ) again. Let one pair of opposite faces carry 
the equations 

Qi(x, xi, xn, x 2 \ , a {2) ) s=1 = 0, Qi{x 3 , x 13 , x 123 , x 23 ; a (2) )s =0 = 0, 
and let two other pairs carry the equations 

Q(x,Xi,x it3 ,x 3 ,a {l \e) = 0, Q(xj, Xij, x 123 , x j)3 , a w , e) = 0, {i, j} = {1,2}, 
where the polynomial 

Q(xi, x 2 , x 3 , x 4 , 7, e) = {xi - x 2 ){x 3 - x 4 ) + ^{e^ 1 - ex 3 x A ) 



actually coincides with (Hf) up to the permutation of x 2 ,x 3 . This awkward structure is 
3D-consistent and, surprisingly, satisfies the tetrahedron property. It can be also extended 
to the lattice Z 3 . 



8 Concluding remarks 

Non-commutative analogues of some of the equations in the Q-list are known. In particular, 
a quantum version of (J^i 5=0) appeared in [27] . In [UJ the consistency approach was 



extended to the non-commutative context, when the fields take values in an arbitrary 
associative algebra. The definition of three-dimensional consistency remains the same in 
this case, however, the assumption on the affine-linearity is replaced by the requirement 
that the equation can be brought to the linear form px = q with respect to any variable x. 
These two properties are not equivalent in the noncommutative case, as is seen from the 
following examples. The first one was found in [UJ and the other two by V.V. Sokolov 
and V. Adler (unpublished): 



a 



a 



(i) 



(x - x 2 )(x 2 - x 12 ) 1 = a {2 \x - x 1 )(x 1 - £12) \ (Oil 5=0) 
(xi — xi 2 + a™)(x — x\ — a*- 1 -*) -1 

= a (2 \x 2 - x 12 + a {2) ){x -x 2 - a (2) )"\ (Qi|«=i) 
(1 - (a (1) ) 2 )(xi - « (2) x 12 )(a (1) x - xx)- 1 

= (1 - (a (2) ) 2 )(a;2 - « (1) x 12 )(a (2) x - x 2 y l . (Q 3 \ s=0 ) 
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The existence of non-commutative analogs of (Q2), ( IQ3I (5=i) and ( Q4 ) remains an open 
question. Although the analysis of the singular solutions may still be useful as a general 
principle, our technique is based on the algebraic properties of affme-linear and biquadratic 
polynomials and is therefore not applicable to this problem. 

More general quantum systems with consistency property were found recently in [HI EJ - 
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